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ABSTRACT 

We present budda (BUlge/Disk Decomposition Analysis), a new code de- 
voted to perform a two-dimensional bulge/disk decomposition directly from the 
images of galaxies. The bulge component is fitted with a generalized Sersic 
profile whereas disks have an exponential profile. No other components are in- 
cluded. Bars and other sub-structures, like lenses, rings, inner bars and inner 
disks, are studied with the residual images obtained through the subtraction of 
bulges and disks from the original images. This means that a detailed structural 
analysis of galaxies may be performed with a small number of parameters, and 
sub-structures may be directly studied with no a priori assumptions. As has been 
already shown by several studies, two-dimensional fitting is much more reliable 
than one-dimensional profile fitting. Moreover, our code has been thoroughly 
tested with artificial data and we demonstrate it to be an accurate tool for deter- 
mining structural parameters of galaxies. We also show that our code is useful in 
various kinds of studies, including galaxies of, e.g., different morphological types, 
and inclinations, which also may be observed at different spatial resolutions. 
Thus, the code has a broader range of potential applications than most of the 
previous codes, which are developed to tackle specific problems. To illustrate its 
usefulness, we present the results obtained with a sample of 51 mostly early-type 
galaxies (but covering the whole Hubble sequence). These results show some of 
the applications in which the code may be used: the determination of param- 
eters for Fundamental Plane and structural studies, quantitative morphological 
classification of galaxies and identification and study of hidden sub-structures. 
We have determined the structural parameters of the galaxies in our sample and 
found many examples of hidden inner disks in ellipticals, secondary bars, nuclear 
rings and dust lanes in lenticulars and spirals, and also wrong morphological 
classification cases. We now make budda generally available to the astronomical 
community. 
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Subject headings: galaxies: fundamental parameters — galaxies: photometry — 
galaxies: statistics — galaxies: structure — techniques: image processing — 
techniques: photometric 



1. Introduction 

It is widely assumed that the light distribution in galaxies may be described by one or 
more components, which might be dynamically and/or chemically distinct. While most of 
ellipticals (E) may be photometrically described by only one major spheroidal component, 
lenticular (SO) and spiral (S) galaxies in general need at least two components: a spheroidal 
bulge and a flattened disk. However, in almost every object we can trace the presence of 
other relevant components, such as primary bars, secondary bars (in spirals and SOs), in- 
ner disks, lenses and rings, none of which have been as carefully studied as the bulges and 
disks. Although there is an unclear and poorly understood interplay between these com- 
ponents, which is not generally taken into account in most structural analyses, bulge/disk 
decomposition has been a major motif behind a number of important contributions to our 
understanding on galaxy formation and evolution. Examples include the Tully-Fisher rela- 
tion (Tully & Fisher 1977), the Fundamental Plane (Faber et al. 1987; Dressier et al. 1987; 
Djorgovski & Davis 1987; Bender et al. 1992, hereafter FP), the correlation between the 
scale lengths of bulges and disks (Courteau et al. 1996; de Jong 1996b; dos Anjos Sz Gadotti 
2004) and the evolution of galaxy morphologies (Marleau & Simard 1998). See also Peng et 
al. (2002) and Graham (2002a,b) for further and more recent examples. 

Two important issues are related to the subject of bulge/disk decomposition: the fitting 
method and the mathematical functions used to describe each component. Generally speak- 
ing, there are two fitting methods, the one-dimensional (ID) and the two-dimensional (2D). 
In the ID method, an azimuthally averaged surface brightness profile of the galaxy under 
study, or the major or minor axis profiles, is fitted by one or more components. This method 
has the advantage of being simple and fast, and works in low S/N conditions. However, it 
has two major problems. Firstly, the procedure to obtain the galaxy profile may vary, and 
each procedure has its own drawbacks. In isophote fitting (Jedrzejewski 1987), which is com- 
monly used, changes in ellipticity and position angle at consecutive isophotes turn the profile 
poorly defined, since it is extracted not along a straight line but an arc. This is especially 
true when a bar or a strong spiral arm is present. Alternatively, one may extract the profile 



^ Based on observations made at the Pico dos Dias Observatory (PDO/LNA - CNPq), Brazil 



-3- 



along the galaxy major axis, which may be itself poorly defined, but in this case much of 
the information in the image is thrown away (Baggett et al. 1998). Moreover, minor axis 
profiles are also essential in this case since the bulge dominates the profile along this axis. 
Secondly, ID fitting averages out in a complicate way non-axisymmetric components, like 
bars. It is possible, though, to obtain reasonable results with ID fitting, as shown by Kent 
& Glaudell (1989), which study the SBO galaxy NGC 936 by means of ID decomposition 
but also incorporating further intricate steps (see also Kent 1986). In 2D fitting, information 
in the whole image is used to build a model for each component and thus differences in the 
projected ellipticities and position angles become an advantage, since they help to constrain 
the parameters of the different components. Moreover, non-axisymmetric components are 
now taken into account. There are several examples in the literature of studies showing that 
the 2D method is much more reliable than the ID (de Jong 1996a), retrieving more accurate 
structural parameters. 

While there are several examples in the literature of studies which use ID fitting (Baggett 
et al. 1998), only in the last decade 2D fitting algorithms become more popular, mainly 
because of the higher computational speed of today's machines. As examples, we may point 
out the algorithms developed by de Jong (1996a), by Simard (1998), and more recently by 
Peng et al. (2002) and Trujillo et al. (2001). Most of the available algorithms were developed 
to tackle specific problems, like, e.g., face-on late-type spirals (de Jong 1996a), distant 
galaxies (Simard 1998), or nearby galaxies observed with very high spatial resolution (Peng 
et al. 2002). On the other hand, it is a common practice to include several components other 
than the spheroid and the disk in an attempt to get better fits. In this way, one can model 
bars and, in the case of Peng et al. (2002), also nuclear point sources, inner disks and inner 
bars, and lenses. In the Peng et al. algorithm 40 free parameters are used, and in several 
cases a single component is fitted by various different functions. 

A caveat behind this philosophy of performing structural analysis is related to the fact 
that we do not know how well a single general function could fit these sub-structures with 
a given accuracy. Therefore these fits might be subject to an unknown systematic effect. 
For instance, as some bars show an exponential profile, especially in late-type galaxies, and 
others have a fiat light distribution, especially in early-type galaxies (Elmegreen 1996), with 
a steep outer drop, it is difficult to foresee how a single general function could be used to 
fit bars and give meaningful structural parameters. This is also the case for lenses and 
other sub-structures. In our code, we explore the alternative route of fitting the two major 
components, representing the disk and the bulge contributions, and to perform a detailed 
study of sub-structures in the residual images obtained after the subtraction of these major 
components. In principle, this procedure could allow us to study the minor structures and 
assess the possibility to describe them as families of functions. It should be noted, however. 
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that residual images can also be obtained with ID fits if one can derive reasonable estimates of 
the apparent axis ratios of the bulge and disk components. For instance, Buta et al. (1992) 
show a residual image of NGC 4622 based on ID decomposition expanded to 2D using 
assumptions about its components projected shapes. Also, 2D bulge/disk decomposition 
is not the only way to bring out hidden features. Using Fourier decomposition (see, e.g., 
Considere & Athanassoula 1988), Buta et al. (1992) were able to bring out remarkably well 
the amazing counter-winding spiral arms in NGC 4622. To a lesser extent, a similar analysis 
can be performed via unsharp masking and color maps (see, e.g., Erwin & Sparke 2002, for 
a recent application). 

Therefore, following these guidelines, we have developed a new and general purpose 
algorithm, budda^ - which stands for BUlge/Disk Decomposition Analysis - devoted to 
perform a detailed 2D structural analysis in galaxy images. The code allows us to perform 
a direct study on sub-structures in galaxies by means of accurate residual images, besides 
determining structural parameters of the main components, namely the bulge and the disk. 
In Sect. 2 we present the details of how the algorithm works. In Sect. 3 we present a 
thorough analysis of the reliability of the algorithm using a set of artificial data images of 
galaxies having known structural parameters. In Sect. 4 we apply our algorithm to a sample 
of 51 galaxies and show the results. Finally, in Sect. 5 we finish with a summary and a brief 
discussion of our results and the potential apphcations of our code. 



From the point of view of our algorithm, a galaxy is simply an image composed by just 
two major components: the bulge and the disk. In order to extract the properties of these 
two components we need to describe them in a generic way and predict their expected con- 
tribution at each pixel of the observed image. For the disk we have adopted the exponential 
profile (Patterson 1940; Freeman 1970): 



where h is the disk scale length, 7o<i represents the central brightness and r is the semi major 
axis of the elliptical isophotes having ellipticity e^. Both components, the bulge and the 
disk, share a common center located at Xo,yo) but the ellipticity (e^, e^) and position angle 



^Scc http://www.astro.iag.usp.br/~dimitri/budda.htnil where BUDDA is available to the astronomical 
community. 



2. The Algorithm 
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{PAd, PAh) are not necessarily equal. In principle, we expect that the information on the 

shape parameters of the disk is mainly sampled in the outermost isophotes of the galaxy, 
while the bulge is more influential in the innermost central region. Therefore, by tracking 
the subtle pixel variations due to changes in the isophotal shape going from the center to 
the external region the program is able to model both the bulge and disk isophotal shapes. 

Bulges are usually described in the literature by the r^^"^ profile or de Vaucouleurs' 
function. Several examples of elliptical galaxies are known to accurately obey the de Vau- 
couleurs' law over as much as 10 intervals of magnitude in surface brightness. On the other 
hand, there arc also several objects, including bulges, where this relation docs not provide 
an accurate description. For these reasons we rather prefer the use of the so called Sersic 
profile, which is a more generic mathematical expression that includes the de Vaucouleurs 
profile as a special case. Moreover, there is strong evidence that the Sersic profile is more 
efficient in the empirical description of isophotal properties of early type galaxies (Caon et 
al. 1993). Therefore, we adopt the solution of describing the brightness profile of the bulge 
by using the relation 

76(r) = /oa0-''"[('-/'-«^)''"l, (2) 

where is the effective radius, is the central brightness and n is an index controlling 
the shape of the brightness profile. For n = 4 we recover the dc Vaucouleurs profile, while 
for n = 1 we have the exponential profile. Therefore the index n is an important quantity 
describing the radial variations of the brightness profile. The numerical constant bn is defined 
in such a way that the the effective radius, reb, contains half of the total bulge luminosity 
and its value is approximately given by 6„ ~ 0.868n + 0.142. At r — Teh, h — hb: the bulge 
effective surface brightness. The coordinate r indicates the semi-major axis of the bulge 
elliptical isophotes centered on XQ,yQ, having an ellipticity ej, and a position angle PAj,. 

The central brightness profile is known to be heavily affected by the atmospheric see- 
ing. Therefore we have included in the code a circular gaussian smearing of the brightness 
profile controlled by a parameter associated with the seeing radius (os). This is similar to 
what was already done by dc Vaucouleurs & Capaccioli (1979) and Capaccioli & de Vau- 
couleurs (1983). Basically it affects only the very central region where the convolution with 
the original profile is able to ffattcn the central brightness distribution. But its inclusion 
is quite important because otherwise we would get very discrepant models in the central 
region, changing the optimum minimum solution. We have also included in the code 
a supplementary correction of the sky level. In fact, all images were sky corrected before 
applying the code. But the fact is that even in a careful analysis the sky subtraction could be 
quite important mainly in the outer region where the disk and bulge brightness profiles come 
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close to zero. An oversubtraction of the sky contribution could easily mimic a truncation 

of the external region of the disk, while an undersubtraction can be misinterpreted as an 
extended disk. Therefore, we have included this correction term {^sky) in order to verify if 
the sky correction was done properly or not. Whenever this correction was higher than a 
few percent of the sky level we have reanalized the sky subtraction procedure to verify the 
reasons behind such discrepancies. In some cases, the sky subtraction might be difficult to 
determine due to the projection of nearby stars, companion galaxies, or simply due to the 
large size of the galaxy in the frame that precludes a good samphng of sky regions. These 
worse cases, where we could not improve the sky subtraction procedure, have been removed 
from our analysis. Note that we have done a careful sky subtraction in the galaxy images 
studied in this paper (see Sect. 4) and that these images of our sample of galaxies generally 
allow for a proper determination of the sky background. Thus, the cases where the A^^^^ 
term suggests that the sky subtraction may have failed were indeed very rare in this study. 

It is also possible to fit edge-on galaxies. In this case, the disk ellipticity is not anymore a 
useful parameter and it is discarded by the code. Instead, it uses a vertical light distribution 
for the disk which follows van der Kruit & Searle (1981), i.e., oc sech^(z/zo), where zq is the 
isothermal scaleheight. Moreover, the ellipses describing the bulge and the disk components 
do not need to be necessarily perfect: boxy and disky isop botes are allowed, and the code 
looks for the best description. These are described by generahzed ellipses and the eUipse 
index parameter 7 controls their shape: for 7 = 1 the ellipses are perfect, while for 7 > 1 
the eUipses are boxy; disky ellipses have 7 < 1. 

The user also has the option of considering an inner truncation in the disk and an 
outer truncation in the bulge. Note, however, that in this study we have not restricted 
the contributions of both bulge and disk in any region of the galaxies whatsoever. We have 
added these truncations as user specifications in the code to allow for the freedom of studying, 
for instance, disks with type II profiles (Freeman 1970). Any user should be very careful, 
though, in using these parameters, by making a priori assumptions that a given area of the 
galaxy is dominated by the bulge or the disk, as was often done with the ID fits. These 
truncations may also work as a way to exclude certain regions from the 2D fit. We stress 
again that this has not been done in this paper. For instance, this is not what we have done 
to deal with galaxies with bars or other conspicuous components. To model such cases, we 
rely on the code's ability to distinguish bars and spiral arms from the bulge and the disk 
by means of differences in the position angle, ellipticity and luminosity profile. Unless these 
three constraints are all equal for, say, the bulge and the bar, the code is able to separate 
these components, as is revealed by the analysis of barred galaxies (see Sect. 4), in which 
the combined luminosity profile of the bulge and disk models falls below the observed profile 
for the galaxy only in the bar region. 
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Altogether the code needs to determine a total number of 11 main parameters: xq, yo, 
hd) h, Cd, PAd, lob, f^eb, Cft, PAh, and n. The others, Ug, Ag^y, the two 7 indices (one for 
bulge and one for disk), and the two truncation radii (the outer one for bulge and the inner 
one for disk) may be set constant in general. An initial guess for these parameters, together 
with an initial allowed range of variation, are provided by the user in an input file read by 
the program. If the range of a given parameter is initialized as zero the code understands 
that it is a fixed parameter that should be kept constant during the iteration process. This 
is a useful way of dealing with parameters that are difficult to be determined even by using 
the whole image information. For instance, if we have an almost circular component the 
position angle becomes badly constrained and difficult to be determined. In these cases 
we have opted to have an initial run of the code, keeping these parameters fixed, and then 
a second run was made leaving them unconstrained. The program also expects to receive 
from the user information about the noise amplitude of the image being analyzed. This is 
a crucial piece of information since the code will decide the best solution by calculating the 
of the overall fit of the model. This information comes in two ways: the fiat noise level 
and the Poisson noise level. The former is a number in counts that represents the noise 
in flat regions of the image, which is dominated by contributions from the noise in the sky 
level and from residuals due to imperfections in the flatfleld normalization. The flat noise 
level can be estimated by the user or by the code itself; the minimum valid count for the 
code will be the square root of the flat noise level. The Poisson noise level is kept to 1 if 
the user wants to assume that the noise in the galaxy image is purely Poissonic. The code 
may also verify if this is a correct assumption and change this level accordingly. To avoid 
complicated regions in the image the user can flag those regions dominated by stars, defects 
or other companion objects that could spoil the solution. This may be done with the IRAF^ 
task IMEDIT. Note that while only bright stars may interfere in the code's determination 
of the bulge and disk parameters, the fainter ones still contribute to the calculation of the 
X^. Thus it is very important that the galaxy image be as free as possible from this kind of 
perturbation. Moreover, the maximum radius of the galaxy can also the informed so that 
the code doesn't try to find information outside the region where the object is present. At 
the end of this initial step the code selects all valid pixels that will be used during the fitting 
procedure. 

The fits are done in raw intensity units and all pixels have the same weight. To obtain 
the optimum solution we have adopted the multidimensional downhill simplex method (Press 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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et al. 1989). The iteration scheme works simply by testing all the possible combinations of 
the parameters in each step and choosing the one having the smallest x^. In each step, 
the multidimensional volume along the n-dimensional parameter space is preserved. If a 
minimum is found along any given direction in the parameter space, the program performs a 
parabolic interpolation to improve the estimate for the next step. When a global minimum 
is found, the program does a final parabolic interpolation and estimates the error along any 
given parameter using the maximum variation allowed by the confidence level demanded 
by the user. Thus, it must be noted that the errors at each parameter determined by 
the algorithm can sometimes be overestimated. This is not related to the fact that the 
decomposition is a 2D one, but only because the code uses the maximum ^ variation allowed. 
In the end, the code generates a file containing the expected image for each component that 
could be added and compared with the original data. 

The code is available in two executable task files: GMODEL and bmodel. The former 
is the one which determines the best bulge and disk models for the galaxy under study; the 
latter builds synthetic images from the results of the former. 

3. Tests 

To verify the reliability of the budda algorithm we have performed a long series of 
tests. Firstly, the code was compiled on different machines (PCs and workstations) and 
Tinder different operational systems (Linux and Unix). As there are small differences between 
Linux and Unix Fortran compilers, we have made small changes in the code so as to have 
proper versions for Linux and Unix. The different versions were tested and it was checked 
that they give the same results in any machine when applied to the same galaxy, whether a 
real object from our sample or a synthetic generated model. 

To certify that the code reliably retrieves the right structural parameters of bulges and 
disks in galaxies, we have applied it to 41 synthetic galaxies. These were built with the 
BMODEL executable within BUDDA so that we know the galaxies' parameters a priori and 
can compare them with the results when we apply the gmodel budda executable to the 
galaxies' images. The synthetic galaxies are made of a bulge and a disk, with no other 
structural component such as a bar, since the main goal of the tests was to verify the 
reliability in the values obtained for the structural parameters of bulges and disks. While 
it is true that the complexity of real galaxies, which contain bars, spiral arms etc., may 
cause difficulties in the retrieving of the correct parameters, we have argued above that 
the code is able to distinguish between the surface brightness contributions from different 
components through constraints such as the ellipticity, position angle and luminosity profile. 
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Thus, bars and spirals are not considered by budda and then appear more conspicuously in 
the residual images. Moreover, when our code is applied on such more complex galaxies, the 
determination of the errors in the estimates of the structural parameters is made accordingly, 
since it takes into account the value of the for the fit (sec Sect. 2). The code was applied 
exactly in the same way as it is done for real galaxies. The first guesses for the structural 
parameters were randomly chosen so that they can initially be quite different from the 
original value used in the simulation. No matter how far from the right value the first guess 
for a parameter may be, the code gradually tends to converge to the right value, within the 
expected error. This shows how robust the results from the code are, even if one doesn't 
have a good estimate about the morphological properties of the galaxy under consideration 
before the algorithm is applied. 

We have done two kinds of tests. In the first one, a good spatial resolution image was 
simulated, similar to most of the studies using ground-based telescopes without active or 
adaptive optics. For these tests we have adopted a formal resolution of a two arcsec seeing 
in a galaxy image with a diameter of around 2 arcmin. The results we present in the next 
section are based on observations generally somewhat better than this. The tests show that 
this resolution is already good enough to determine structural parameters with excellent 
reliability (see below). To check how BUDDA works in worse cases (e.g., imaging of dwarf 
galaxies or galaxies at intermediate redshifts), we applied the same tests with a resolution 
two times poorer than that. 

It is worth noting that in these tests there was no user intervention, as should be the 
case with real galaxies, when some actions are recommended to be taken in advance in order 
to maximize the code efficiency in retrieving accurate data. For instance, in real galaxies 
GMODEL should be apphed several times, choosing different sets of initial values, to search for 
the best result. That precaution avoids the adoption of a solution that could be located 
in a local minimum. In some cases, one may initially fix some geometric parameters by other 
means; for instance, we could guess the initial position angle or ellipticity with the ellipse 
task from IRAF. Then, after a first solution is found, one may leave these parameters free to 
vary in a second run, using the initial values as those of the solution of the first run. This is 
the recommended procedure, for instance, when the ellipticity is too small (around 0.1) and 
the position angle is difficult to determine. However, none of these actions was done in the 
tests, so that their results are actually a lower limit for the code's accuracy. Poissonic noise 
was added in some of the synthetic images but we found it to be a minor source of error, as 
long as the S/N ratio is reasonably high (i.e., around 50 in the central regions of the galaxy). 

The results from the tests are shown in Fig(s). 1 and 2 for the good resolution case, 
and Fig(s). 3 and 4 for poor resolution. From Fig. 1 we can observe that the reliability 
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of the code is excellent in recovering the surface brightness, scale length and position angle 
both of the bulges, indicated by the open squares, and disks, indicated by the open circles. 
While there are some minor trends, the agreement is always within the estimated fitting 
error delivered by the program. The ellipticities of bulges and disks are more difficult to 
estimate and are subject to larger errors. The reason is that for a fixed value of ellipticity, 
within the formal error, we can always find a set of the other photometric parameters that 
will result in a final error close to the best model. Therefore, the final error is relatively 
insensitive to an error in the intrinsic ellipticity of the two components. Moreover, this effect 
tend to be worse when the bulge/disk ratio is larger than 1, since in that case a change in the 
disk eUipticity practically has no impact in the total error budget of the model. Therefore, 
great care should be exercised when interpreting the ellipticity results. The ellipticity of the 
most eccentric bulges may be underestimated, but again is within the errors. This ellipticity 
issue may be easily corrected with the actions discussed in the above paragraph, by using, 
for instance, IRAF tasks to have an initial value based on the observed image. The solid 
symbols in Fig(s). Id and 3d show how the code works in the cases when the ellipticities of 
bulge and disk are treated as in real galaxies. The Pearson correlation coefficient R in these 
figures are relative only to the sohd symbols. In Fig. 2 we present the results regarding the 
Sersic index and the bulge/disk ratios. For large values of n the adopted seeing can cause 
an important variation in the estimated value. However, for the vast majority of galaxies, 
including Es, this effect is not important since that index is bound to have a lower value 
close to the de Vaucouleurs {n = 4) model. The good agreement also found in Fig(s). 3 and 
4, where we have tested the impact of having low spatial sampling, is encouraging. 

To verify the error produced in the residual images by the differences in the input and 
retrieved synthetic models, we have subtracted the latter from the former in Fig. 5. The 
left panels show of good spatial resolution while the right ones show of poor 

resolution. They show that the discrepancy is restricted to the central region, which has a 
size of the order of the seeing. The right panels represent an especially bad fit when the 
ellipticities were not treated as must be done in real galaxies. Such a problem may be easily 
identified by the bipolar pattern left. The entire original synthetic galaxy occupies the whole 
frame in the upper panels. The brightness and contrast levels displayed were chosen as to 
enhance the differences. 

Another interesting point regarding these tests is that the code is able to identify disks 
which are up to 6 times fainter than the bulge (see Fig(s). 2 and 4), although this value 
depends somewhat on the geometric differences between disk and bulge (in fact, our results 
below show that the code can be even more sensitive). In any case, this result reinforces the 
BUDDA finding of Gadotti & de Souza (2003) on the absence of disks in the SBO galaxies NGC 
4608 and 5701, since the bulge/disk ratio in normal SOs is around 2 (Binney & Merrifield 
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1998). 

As noted above, the errors determined by the algorithm may be sometimes overesti- 
mated. This seems to be especially true for the brightness and for the position angle. The 
1(7 error bars for the surface brightness are typically around 0.75 mag. However, from Fig(s). 

1 and 3, one sees that fiducial errors may be considered to be around 0.25 mag. In these 
figures, there are no error bars associated with the position angles because most of the time 
the algorithm is not able to estimate them reliably, a fact that is related to the ellipticity 
issue discussed above. There are no error bars also related to the bulge/disk ratios (Fig(s). 

2 and 4), but this is because they are calculated via the model images. 

4. A Structural Analysis Study of 51 GalzLxies 

In order to test budda in real galaxies we have observed 51 galaxies at the Pico dos 
Dias Observatory (PDO/LNA - CNPq, Brazil). The only criterion these galaxies obey is 
to be brighter than around 14 in B, and emphasis was put on Es and SOs. The CCD 
observations were done with a 24 inch telescope having a focal ratio f/13.5, and using a 
thin back-illuminated CCD SITe SI003AB, with 1024 x 1024 pixels. The plate scale is 0.57 
arcsec/pixel, resulting in a field of view of approximately 10 x 10 arcmin. The CCD gain 
was set at 5 e^'/ADU and the read-out noise on 5.5 e~. All objects were observed in the R 
passband of the Cousins system. For each object, we have done 3 exposures of 300 seconds. 
The multiple exposures aim to ease cosmic ray removal. The data were calibrated with a set 
of standard stars from (Graham 1982) and corrected for atmosphere and Galactic extinction. 
The later correction was done using the maps of Schlegel, Finkbeiner & Davis (1998). The 
standard processing of the CCD data includes bias subtraction, fiatfielding and cosmetics. 
The first step in the sky subtraction was done by editing the combined images, removing the 
galaxy and stars. After that step we determined the mean sky background and its standard 
deviation (cr). Then, we removed all pixels whose values were discrepant by more than 3 
a from the mean background. A sky model was obtained by fitting a linear surface to the 
image, and this model was subtracted from the combined image. We finally removed objects 
such as stars and Hll regions. All these procedures were done using the IRAF package. The 
observations were done in 6 nights in October 1997. Typical values for the seeing and zero 
point error are, respectively, 0.8 to 1.2 arcsec and 0.02 to 0.05 mag. 

Besides testing BUDDA, these data also allow us to perform a detailed structural analysis 
of galaxies. Since our primary concern is to verify sub-structures in early-type galaxies, our 
sample contains mostly Es and SOs, becoming also suitable to Fundamental Plane studies. 
However, our sample spans the whole Hubble sequence, with some examples of Ss, which are 



- 12 - 



useful to test budda in galaxies with spirals, bars, dust lanes and morphological peculiarities. 
Figure 6 shows the results for all galaxies. For each of them we have obtained with budda the 
structural parameters and artificial images of the model galaxy, as well as of the model bulge 
and disk, separately. Using the IRAF ELLIPSE task, we built the surface brightness profiles of 
the galaxy, of the model galaxy, and of the bulge and the disk, as well as ellipticity, position 
angle and the b4 Fourier coefficient profiles for the galaxy and for the model galaxy. The 
behavior of the b4 coefficient indicates how much the isophotes deviate from a perfect eUipse. 
A positive value for b4 indicates isophotes with a disky component, whereas a negative value 
indicates boxy isophotes. If b4 = 0, then the isophote is a perfect eUipse. This whole 
procedure allows us to evaluate the quality of the fit and other structural parameters of 
the galaxy. Moreover, we obtained residual images, subtracting from the galaxy the model 
images, to verify the presence of sub-structures. In Fig. 6 we show, from left to right, and 
from top to bottom: 

• Original Images, passed through a median filter and converted to surface brightness, 
with the isophotal map (with a 0.5 mag interval between each contour level) 

• Total Residual Images, i.e.. Original Image divided by Complete Model^ 

• Disk Residual Images, i.e.. Bulge Model subtracted from Original Image 

• Bulge Residual Images, i.e.. Disk Model subtracted from Original Image 

• Surface Brightness profiles in the R-broadband, in magnitudes per squared arcsecond, 
as a function of the radius, in pixels, where 

— points with error bars represent the Original Image 

- the short-dashed fine represents the Disk Model 

- the long-dashed fine represents the Bulge Model 

- the full line represents the Complete Model 

— points with full line represent the Total Residual profile (original image — com- 
plete model -|- 28), in magnitudes per squared arcsecond 

• ellipticity, position angle, and b4 Fourier coefficient profiles, where 

— points with error bars give values from Original Images 



''Note that the total residual image is a ratio image, in which discrepancies are emphasized. 
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- the full line represents the Complete Model 

Through a careful analysis of Fig. 6, one can identify and study several features which 
would otherwise remain hidden in the direct galaxy image. Some points must be stressed 
concerning the analysis of residual images. In the computer, one may always change the 
brightness and contrast levels to carefully identify what is real and what is spurious. To 
make the figures, these levels are chosen to highlight features one wants to point out. Due to 
seeing effects the first few pixels in the nucleus are always brighter than the model and this 
excess extends to a few more pixels because there are too few data points for the model to 
be accurate. Thus, analysis in the very inner parts of the galaxy is always dubious. A level 
of subjectivity still remains in this kind of analysis, and because of that, the existence of any 
sub-structure may be confirmed also through the analysis of the radial profiles. For instance, 
an inner disk or bar may reveal itself in surface brightness radial profiles as a bump. Some 
other signatures are described in the following subsections. 

Regarding our results, a first point that should be mentioned is that the broadest clas- 
sification one may do with our sample is to distinguish between galaxies with only one 
spheroidal component, and galaxies that contain also an exponential disk. Behind this ap- 
parently banal statement, there are elliptical galaxies, which shouldhe pure spheroidals, but 
contain inner disks, and lenticular galaxies, which should have a disk, but do not show any 
discernible disk, either because of wrong classification or secular evolution (in the case of 1 
SBO, namely NGC 2217). Moreover, we find several other cases of wrong morphological clas- 
sification and many galaxies with hidden sub-structures. We discuss these generic cases in 
the following subsections. Table 1 presents our results, including the structural parameters 
found for each galaxy. 

We have made an extensive search in the hterature for similar estimates for the galaxies 
in our sample, but it turned out that only to one galaxy a 2D decomposition code has been 
applied before in the R broadband. This is NGC 7177, studied also by de Jong (1996a). 
Unfortunately, however, our bulge parameters can not be compared since de Jong a priori 

assumed an exponential bulge for this galaxy, while our results show that its bulge is better 
described by a Scrsic profile with n = 3.52. On the other hand, the disk parameters may be 
compared. The agreement shown is good indeed: while de Jong's results point to Iq^ = 18.73 
R mag arcsec"^ and h = 15.99 arcsec, BUDDA estimates loa = 18.58 R mag arcscc"^ and 
h = 12.97 arcsec. It is interesting that, at least in this case, using a different bulge luminosity 
profile did not alter significantly the estimate of the disk parameters. Moreover, the fact 
that de Jong also used a model for the bar of this galaxy indicates that BUDDA was able to 
distinguish the contribution to the total luminosity from the bulge, disk and bar, without 
assuming any luminosity distribution for the bar, which is very encouraging. 
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4.1. Inner Disks in Ellipticals 

There are 23 galaxies in our sample classified as Es in the RC3 (de Vaucouleurs et al. 
1991). From these, only 9 seem to us to be pure spheroids, since the fit is excellent with only 
the spheroidal component and the residual images do not show any sub-structure (except in 
the nucleus, where seeing prevents any analysis). Prototypical examples are NGC 822, 6578, 
6958, 7507 and 7619. Others that fall in this category are: IC 4842, 4943, NGC 720, and 
1199. But even on some of these examples some doubt may be cast regarding the presence 
of any more complex morphology. Burstein et al. (2001) suggest that a// galaxies have disks 
but 20% of them are not forming stars and a small fraction are too faint. 

Eight galaxies classified as Es do show inner discs. These may be identified directly by 
the code, which in this case finds the exponential profile. This is the case for e.g. NGC 
1052 and 7192. In some cases, however, the disk is too faint or too small to be recognized by 
BUDDA but it appears clearly in the residual images. One must bear in mind that these inner 
disks are not necessarily exponential disks, which is what our code tries to fit. Examples 
of this case are ESO 115-G008, PGC 64863 and NGC 821. Other signatures of the inner 
disk are disky (b4 > 0) inner isophotes and a distinct behavior of the ellipticity and position 
angle radial profiles from the inner to the outer regions of the galaxy. Disky isophotes are 
a signature of inner disks because it was found that Es with disky isophotes rotate faster 
(Kormendy & Djorgovski 1989). However, the isophote being disky only means that there is 
an excess of light along the major and minor axes of the galaxy, which makes it look pointy. 
Bars seen close to or face-on also may produce disky isophotes. When b4 < 0, the isophotes 
are boxy. This means that the galaxy has a deficiency of light along its main axes, which 
causes a rectangular shape. Boxy isophotes are also found in some Es, which have negligible 
rotation (Kormendy & Djorgovski 1989). Finally, boxy isophotes are also related to bars, 
but when seen edge-on (Binney & Merrifield 1998). Other Es with inner disks are: NGC 
1209, 1700 and 7196. 

The remaining 6 galaxies are examples of bad morphological classification and will be 
discussed below. In all these galaxies we have found a large exponential disk which puts 
them in the lenticular class. In the misclassification section there are also galaxies classified 
as SOs, but that are actually Es, some of them with inner disks or complex morphology. 

4.2. Sub-structure in Lenticulcir Galetxies and another SBO without Disk 

There are 20 galaxies in our sample classified as SOs in the RC3. Eleven are genuine 
lenticulars, in which budda determined the disk parameters, except in one, namely NGC 
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2217. This galaxy apparently has a very strong bar that could have destroyed the pre-existent 
disk within the bar length, and may be another case as those (NGC 4608 and NGC 5701) 
discovered by Gadotti & de Souza (2003); see also Athanassoula (2003). Note, however, that, 
while NGC 4608 and NGC 5701 have indeed strong bars, as denote by their Fourier m = 2 
component (Gadotti & de Souza 2003), it remains to be seen that the bar in NGC 2217 is 
really strong. As pointed out by Buta & Block (2001), bar strength can be defined in terms 
of tangential forces relative to mean radial forces, and the latter are certainly significant 
in galaxies with large bulges, such as NGC 2217. Note, furthermore, that the existence of 
conspicuous rings outside the bar in these galaxies may be an indication of the existence of 
an outer disk, as the bar is only able to consume the disk within it. 

Three of the eleven genuine lenticular galaxies are strongly barred and there are also 
three more strongly barred SO/a galaxies in our sample. In many of these systems one 
encounters interesting sub-structures which do not appear in their direct images. In some 
cases, structural peculiarities may be seen and recognized in a careful analysis of the direct 
images, but are certainly better defined in the residual ones. In NGC 1326 and 2217, the 
inner ring surrounding the bar, already identified by de Vaucouleurs & de Vaucouleurs (1964), 
appears more clearly in the residual image (for a detailed analysis of NGC 1326, see Buta 
et al. 1998). Surprising dust lanes are visible in NGC 1316 and 1947, while the residual 
image of NGC 7280 displays a strong lens. An extremely interesting case is NGC 128. This 
edge-on galaxy, which is a prototype of peanut-shaped bulges was fitted by our code using 
its edge-on capabihties (see Sect. 2). The other objects which fall in this category are: IC 
5250A, NGC 467, 2293, 6673, and 7377. The bulge/disk luminosity ratios for the lenticulars 
in our sample are around 2 with some noticeable deviation cases, in agreement with Binney 
& Merrifield (1998). 

In the remaining 9 galaxies we haven't found a large disk. These are clearly Es rather 
than SOs and will be discussed as misclassification cases below. 



4.3. Wrong Morphological Classification 

To distinguish between an elliptical and a lenticular galaxy is a hard task with visual 
inspection only, especially for faint and/or distant galaxies. Thus confusion between Es and 
SOs is common and may badly interfere in studies of e.g. scaling laws and the origin of the 
Hubble sequence. This is a point where codes like BUDDA may play an important role given 
their higher reliability. In our sample, 6 Es were found to have large disks which means they 
are actually SOs: IC 1459, NGC 1172, 1549, 2305, 7145, and 7778. In 3 of them (IC 1459, 
NGC 1549 and 7145) a lens is clearly visible in the residual images. 
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On the other hand, 9 galaxies classified as SOs were found to have no large disk, which 
puts them in the E bin. As was also the case for some genuine ellipticals, some of these 9 
systems show a morphology more complex than a pure spheroid. Prominent cases are the 
disturbed morphology of NGC 7252, the mner disk in NGC 2271, and the lens in NGC 1553, 
whose ratio image displays also an inner spiral in the center. The last case is especially 
intriguing since lenses are common in SOs. Considering the results from our code, either this 
is a rare case of a lens in an elliptical galaxy, or (more likely) budda was not able to find 
the disk in this lenticular galaxy simply because its disk has not the common properties of 
exponential disks generally found in other disk galaxies. Its disk could possibly be inner 
truncated or have a fiat inner profile (type II profile). The fact that its lens is also very 
bright and has an unusual appearance may have prevented our code to naturally find reliable 
estimates for the disk structural parameters. Indeed, BUDDA finds a disk in this galaxy, but 
that is a very faint disk with unusually high error estimates. It is worth noting that, as 
the position angle and the ellipticity of the lens are very similar to those of the bulge (and 
possibly the disk), it is very hard to disentangle the contribution of the lens to the total 
luminosity of this galaxy without a priori assumptions. NGC 1553 is certainly a case in 
which a more detailed analysis is in order. The other misclassified Es are: IC 1919, NGC 
541, 2205, 6849, 7289, and PGC 67633. 

It is interesting to point out that from the 43 Es and SOs galaxies in our sample (as 
considered in the RC3) 15 were shown to be in the wrong morphological bin, i.e., ^ 35%. 
Moreover, it seems (not surprisingly) that it is equally easy to misclassify an elliptical as an 
SO and the other way around. 

4.4. Sub-structure in Spiral Galeixies 

Eight galaxies in our sample have a morphological class later than SO according to the 
RC3. In many of them hidden sub-structures emerge from the residual images. We point 
out: the bar and the strong spiral arms in NGC 1079; the amazing case of NGC 1291 with 
the nuclear bar and inner ring^ that only show out clearly in the residual images [but note 
that de Vaucouleurs (1975) already recognized these features]; the bar and inner ring of NGC 
7177; the bar and the unusual spiral arms of NGC 7371, and finally the weird spiral structure 
of NGC 7824. This last case is especially puzzling because, in spite of the spiral structure, a 



^This ring may be actually a lens. Since a lens is like a ring of low surface brightness amplitude (see 
Buta & Combes 1996), when a declining bulge and disk background is removed the lens then becomes very 
ringlike. In this case this may be particularly true given the broad appearance of this structure (but note 
that it is connected with the ansae at the ends of the bar). 
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model containing only a bulge is able to fit the surface brightness profile very well. Again, 
this may be a case where the disk can not be modelled as it is traditionally done, i.e., with 
a exponential disk. The inner spirals seem to be detached from the outer ones. And the 
latter have an unusual broad appearance. Bars also reveal themselves in radial profiles of 
ellipticity (monotonically growing) and position angle (constant). The other spirals in our 
sample are: NGC 1637, 7070A, and 7171. 

It is worth noting that in cases like NGC 1079 and NGC 1291, where our images do 
not sample the outermost parts of the galaxy, a better result for the disk model might be 
achieved with deeper and larger images (especially because, in these examples, the fits are 
affected by the outer ring feature). Also, in NGC 7171, the bulge looks overly important 
at large radii, which may be causing the bright halo around the ratio image. This is a 
case in which a more realistic model may be accomplished at the expenses of some a priori 
assumptions, like e.g., keeping the Sersic index constant or assuming an outer truncation for 
the bulge. 



4.5. Nucleetr Activity and Companions 

There are 7 galaxies in our sample which are classified as AGNs at some level (e.g.. Liner, 
Seyfert) in the Nasa Extragalactic Database^. They are: IC 1459, NGC 1052, 1316, 1326, 
1553, 7177, and 7280. There are also 10 galaxies which are classified as multiple systems 
in the Leda Extragalactic Database^ or show a clear companion in our images. This does 
not mean galaxies in groups or clusters and does not exclude systems that are only due to 
projection. It only means that these galaxies have a close companion in the plane of the sky. 
These are: ESQ 115-G008, IC 5250A, NGC 128, 467, 1052, 2293, 7070A, 7196, 7252, and 
PGC 67633. The only galaxy which is a member of both sub-samples is NGC 1052. We have 
not found any clear tendency in these 2 groups of galaxies that makes them distinct from 
the whole sample. Only two wrong classification cases have companions, which means that 
these cases may not be caused by some unknown tidal process that could slightly change the 
morphology of a galaxy. 



^Seehttp: / /nedwww .ipac.caltech.edu/ 
''See http://leda.univ-lyonl.fr/ 
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5. Concluding Remeirks 

We have thoroughly presented BUDDA, a code devoted for bulge/disk decomposition in 
galaxy images. As was shown in details above, its development and performance tests show 
that it is a very suitable tool for the structural analysis of galaxies. The parameters for the 
bulge and disk models can be retrieved with great confidence as well as the residual images 
obtained from them. Its simplicity to use and its broad range of applications add to its 
qualities. 

Applying BUDDA on our sample of 51 galaxies we could check how it works concerning 
its many potential applications. One obvious use of the code is to obtain bulge and disk 
structural parameters to study e.g. the Fundamental Plane and other scaling relations. 
Such a study regarding the galaxies in our sample will appear in a forthcoming paper. Such 
relations may also be strengthened if wrong morphological classification cases are identified, 
which is reliably done with our code, as we showed in Sect. 4.3. Galaxies classified as Es, 
but actually being SOs, which was the case for 6 of our galaxies, may add to the spread 
in scaling relations. Obviously, the opposite is also true. In the FP, for instance, SOs have 
generally lower mass-to-light ratios than Es (Bender et al. 1992). 

Another use of BUDDA, in which we've found it to be very powerful, is to reveal hidden 
sub-structures, especially in the inner regions of galaxies. In Sect(s). 4.1 and 4.3 we revealed 
inner disks in many elliptical galaxies, which may be the result of a recent merger. Or are 
they related to galaxy formation processes? BUDDA has found also several examples of inner 
spirals, inner and nuclear rings, as well as inner and secondary bars in SOs and Ss. In some 
cases a hint of them shows up in the direct images or in the profiles, but in others there is no 
a priori sign of them at all. These sub-structures, which appear conspicuously in the residual 
images, are relevant for e.g. dynamical studies, since they are related to non-circular orbits 
and Lindblad resonances. Moreover, these sub-structures may possibly be relevant to the 
long standing problem of fueling active galactic nuclei (Combes 2001), and may be directly 
studied in the residual images. 

We thank the Pico dos Dias Observatory staff for support during the observing runs. 
Financial support from FAPESP grants 99/07492-7 and 00/06695-0 is acknowledged. We 
would like to thank the anonymous referee for many enlightening comments and helpful 
suggestions. 



-19- 



REFERENCES 

Athanassoula, E. 2003, MNRAS, 341, 1179 

Baggett, W. E., Baggett, S. M., & Anderson, K. S. J. 1998, ApJ, 116, 1626 
Bender, R., Burstein, D., & Faber, S. 1992, ApJ, 399, 462 
Binney, J., & Merrifield, M. 1998, Galactic Astronomy, CUP 

Burstein, D., Saglia, R. R, CoUess, M., Davies, R. L., McMahan, R. K., & Wegner, G. 2001, 
in Galaxy Disks and Disk Galaxies, ASP Conf. Series, 230, 153 

Buta, R., Alpert, A. J., Cobb, M. L., Crocker, D. A., & Purcell, G. B. 1998, AJ, 116, 1142 

Buta, R., & Block, D. L. 2001, ApJ, 550, 243 

Buta, R., & Combes, F. 1996, Fund. Cosmic Phys., 17, 95 

Buta, R., Crocker, D. A., & Byrd, G. G. 1992, AJ, 103, 1526 

Caon, N., Capaccioli, M., & D'Onofrio, M. 1993, MNRAS, 265, 1013 

Capaccioli, M., & de Vaucouleurs, G. 1983, ApJS, 52, 465 

Combes, F. 2001, in Advanced Lectures on the Starburst-AGN Connection, Proceedings of 
a conference held in Tonantzintla, Puebla, Mexico, 26-30 June, 2000. Edited by Itziar 
Aretxaga, Daniel Kunth, and Raul Mujica. Singapore: World Scientific, 2001., p. 223 

Considere, S., & Athanassoula, E. 1988, A&A, 76, 365 

Courteau, S., de Jong, R., & Broeils, A. 1996, ApJ, 457, L73 

de Jong, R. S. 1996a, A&A, 118, 557 

de Jong, R. S. 1996b, A&A, 313, 45 

de Vaucouleurs, G. 1975, ApJS, 29, 193 

de Vaucouleurs, G., & Capaccioh, M. 1979, ApJS, 40, 699 

de Vaucouleurs, G., & de Vaucouleurs, A. 1964, in: Reference Catalog of Bright Galaxies, 
Austin: University of Texas 



-20- 



de Vaucouleurs, G., de Vaucouleurs, A., Corwin, H. G., Buta, R. J., Paturel, G., & Fouque 
P. 1991, in: Third Reference Catalog of Bright Galaxies, Springer- Verlag, New York 
(RC3) 

Djorgovski, S, & Davis, M. 1987, ApJ, 313, 59 
dos Anjos, S., & Gadotti, D. A. 2004, submitted 
Dressier, A., et al. 1987, ApJ, 313, 42 

Elmegreen, D. M. 1996, in Barred Galaxies, ASP Conf. Ser., 91 (San Francisco: ASP); edited 
by R. Buta, D. A. Crocker & B. G. Elmegreen, p. 23 

Erwin, P., & Sparke, L. S. 2002, AJ, 124, 65 

Faber, S., et al. 1987, in Nearly Normal Galaxies, ed. S. M. Faber (New York: Springer), 
175 

Freeman, K. C. 1970, ApJ, 160, 811 

Gadotti, D. A., & de Souza, R. E. 2003, ApJ, 583, 75 

Graham, A. W. 2002, MNRAS, 334, 859 

Graham, A. W. 2002, MNRAS, 334, 721 

Graham, J. A. 1982, PASP, 94, 244 

Jedrzejewski, R. I. 1987, MNRAS, 226, 747 

Kent, S. M. 1986, AJ, 91, 1301 

Kent, S. M., & Glaudell, G. 1989, AJ, 98, 1588 

Kormendy, J., & Djorgovski, S. 1989, ARA&A, 27, 235 

Marleau, F. R., & Simard, L. 1998, ApJ, 507, 585 

Patterson, F. S. 1940, Harvard Obs. Bull, 914, 9 

Peng, C. Y., Ho, L. C, Impey, C. D., & Rix, H-W. 2002, AJ, 124, 266 

Press, W. H., Flannery, B. P., Teukolsky, S. A. & Vetterling, W. T. 1989, in Numerical 
Recipes, CUP 

Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 



- 21 - 



Simard, L. 1998, in Astronomical Data Analysis Software and Systems VII, A.S.P. Conference 
Series, Vol. 145, 1998, R. Albrecht, R.N.Hook and H.A. Bushouse, eds., p.l08 

Trujillo, L, Aguerri, J. A. L., Gutierrez, C. M., & Cepa, J. 2001, A J, 122, 38 

TuUy, R. B., & Fisher, J. R. 1977, A&A, 54, 661 

van der Kruit, R C, & Searle, L. 1981, A&A, 95, 105 



This preprint was prepared with the AAS IM^jX macros v5.0. 



-22- 



Fig. 1. — Comparison between the parameters of the 41 synthetic galaxies as retrieved by 
BUDDA in the case of good spatial sampling: (a) central surface brightness of the disk and 
effective surface brightness of the bulge in arbitrary magnitudes; (b) characteristic radius 
of the disk and effective radius of the bulge in arbitrary units; (c) position angle, and (d) 
ellipticity. Although in (a) and (b) the units are arbitrarily chosen, all values displayed here 
are similar to those found for nearby observed galaxies. Circles are for the disc and squares 
are for the bulge parameters. The error bars were not considered for clarity reasons but 
a typical error bar is shown in each panel. Note that the typical error bar considered is 
as determined by BUDDA, which is not necessarily similar to the scatter in this figure. For 
instance, in (a) it is possible to see that our code typically overestimates the errors for the 
intensities. In (c) there are no error bars, for the code typically does not find a meaningful 
error estimate for the position angles, since variations in these parameters do not cause 
significant changes in . The Pearson correlation coefficient R of each test is also indicated. 
For (d), R is relative only to the solid symbols. These refer to those tests in which the first 
guesses for the ellipticities were not random, but determined from the results of the ellipse 
task in IRAF, as is done for real galaxies. 

Fig. 2. — Same as Fig. 1 but for (a) Sersic index, and (b) the bulge/disk ratios in luminosity 
(circles) and size (squares). 

Fig. 3. — Same as Fig. 1 but for poor spatial samphng. Typical error bars are as in the 
good samphng case (see Fig. 1). 

Fig. 4. — Same as Fig. 2 but for poor spatial sampling. Typical error bars are as in the 
good sampling case (see Fig. 2). 

Fig. 5. — Examples of the discrepancy between a synthetic galaxy built with BMODEL and 
the one retrieved by GMODEL within BUDDA. The upper left panel shows a case of good 
spatial resolution while the upper right one shows a case of poor resolution. The last is an 
example of a specially bad fit. Even so, the differences are restricted to the central region, 
since the entire original galaxy occupies the whole frame. The brightness and contrast levels 
displayed were chosen as to enhance the differences. The lower panels show in details the 
central region of each of the former images, respectively. 

Fig. 6. — Results of the BUDDA algorithm applied in R-band CCD images for all the 51 
galaxies in our sample. At the top it is shown the name of the galaxy, its morphological 
type according to the RC3 and the size of the images in arcseconds. See the text for further 
explanations. The remaining figures are available in the electronic edition of the Journal. 
The printed edition contains only a sample. 
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Table 1. Results from the structural analysis for 51 galaxies. 
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NGC2217 
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22.24 


59.01±7.15 
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no disk 


NGC2271 


SABO 










22.00 


28.94ibl.02 


2.81±0.25 


trueE+inner disk 


NGC2293 


SABO(s) 


1.41 


1.91 


21.87 


72.14±17.57 


21.47 


28.67±2.46 


4.06±0.94 




NGC2305 


E2 


2.80 


1.84 


22.49 
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21.03 


19.75±1.09 


2.99±0.45 
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NGC6578 


E2 










20.80 


43.11±2.58 


4.84±0.72 




NGC6673 


SABO(s) 


7.92 


2.06 


20.69 


33.72ib34.14 


20.19 


34.60±3.16 


5.07±1.11 




NGC6849 


SBO 










21.73 


51.57±2.62 


4.12it0.42 


trueE 


NGC6958 


E2 










20.18 


36.98±2.17 


4.49±0.68 




NGC7070A 


Irr 


1.26 


1.64 


19.75 


20.43d=2.73 


20.92 


22.65±2.04 


2.74±0.48 


dust lane 


NGC7145 


EO 


4.65 


1.44 


20.88 


41.37d=7.14 


21.06 


44.35±1.25 


4.88±0.29 


trueSO+lens 



Table 1 — Continued 



Galaxy 


T 


B/D 


B/D 






lab 






Comments 






(Luiinii(jsi!\') 


(Size) 


[H liiag aicscc " 


I (arcscr) 


[R mag aicscc " 


j (aix-scc) 






NGC7171 


SBb(rs) 


0.44 


1.12 


20.18 


41.28±8.00 


21.25 


40.57±10.18 


2.17±1.38 




NGC7177 


SABb(r) 


7.55 


1.63 


18.58 


12.97±4.07 


20.19 


34.89it2.86 


3.52±0.68 


inner ring 


NGC7192 


EO 


118.52 


1.44 


21.16 


7.98±10.44 


20.61 


39.05±1.22 


4.06±0.30 


inner disk 


NGC7196 


E2 










20.05 


30.65±1.45 


4.05±0.48 


inner disk 


NGC7252 


SO 










19.92 


19.72±0.99 


4.03±0.55 


trueE? 


NGC7280 


SABO(r) 


2.96 


2.07 


20.34 


27.47±4.94 


20.76 


28.78±1.41 


5.32±0.67 


lens 


NGC7289 


SO(r) 










21.86 


33.61ib2.15 


4.56±0.63 


trueE 


NGC7371 


SO/a(r) 


1.10 


1.46 


20.00 


35.49±2.48 


21.98 


35.72±1.72 


3.27±0.37 


bar 


NGC7377 


SO(s) 


2.62 


1.63 


21.09 


75.53±19.66 


20.87 


52.63±2.00 


4.01±0.33 




NGC7507 


EO 










20.08 


46.97±1.46 


4.35±0.35 




NGC7619 


E2 










20.51 


41.79±0.78 


4.34±0.19 




NGC7778 


EO 


1.34 


1.36 


21.10 


55.86±0.29 


19.75 


11.40±0.29 


3.40±0.10 


trueSO 


NGC7824 


Sab 










20.70 


21.32±1.02 


4.47±0.50 


no disk? 


PGC64863 


E4 










21.33 


19.37±1.43 


3.29±0.50 


inner disk 


PGC67633 


SO 










21.80 


21.88±2.05 


4.32±0.88 


trueE 



"^Errors shown correspond to la. 



